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ABSTRACT 

Context. Understanding the depletion of heavy elements is a fundamental step towards determining the structure of pre- 
protostellar cores just prior to collapse. 

Aims. We study the dependence of the NO abundance on position in the pre-protostellar cores L1544 and L183. 
Methods. We observed the 150 GHz and 250 GHz transitions of NO and the 93 GHz transitions of N2H+ towards L1544 and 
L183 using the IRAM 30 m telescope. We compare the variation of the NO column density with position in these objects with 
the H column density derived from dust emission measurements. 

Results. We find that NO behaves differently from N2H"'" and appears to be partially depleted in the high density core of L1544. 
Other oxygen-containing compounds are also likely to be partially depleted in dense-core nuclei. The principal conclusions are 
that: the prestellar core L1544 is likely to be 'carbon-rich'; the nitrogen chemistry did not reach equilibrium prior to gravitational 
collapse, and nitrogen is initially (at densities of the order of 10* cm~'^) mainly in atomic form; the grain sticking 
probabilities of atomic C, N and, probably, O are significantly smaller than unity. 
Conclusions. 
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1. Introduction 

It is now accepted that most, if not all, C-containing 
species deplete on to dust grain surfaces in the central, 
high-density regions of prestellar cores. Thus, when one 
compares a map of C^^O, for example, with that of the 
dust emission, one sees that the CO isotopomer traces 
only an outer shell and not the high-density interior 
(riH <; 3 X 10^ cm~^)^; this limits our understanding of 
prestellar cores, because the molecular line emission pro- 
vides information on the kinematics. 

Surprisingly, N-containing species (in particular, NH3 
and N2H+) manage to remain in the gas phase up 
to densities of a few times 10^ cm^"^. It is not en- 
tirely clear what happens at still higher densities, but 
it seems likely that the N-containing species disap- 
pear from the gas phase when nn ^ 10^ cm^'^ 



Send offprint requests to: C.M. Walmsley 

^ In this paper, we express gas densities implicitly in terms 
of the total number density of protons, uh = ?i(H) + 2n(H2) + 
n(H"'") + ... ~ 2n(H2). However, to be consistent with some of 
the cited publications, we quote explicitly values of n(H2) from 
time to time. 



(see, for example, Belloche & Andre 2004, Pagani et 
al. 2005). The different behaviour of the N-containing 
species was originally attributed to the relatively high 
volatility of N2 on grain surfaces ( [Bergin fc Langer 1997 
Aikawa et al. 2001 ). but recent laboratory measurements 



dOberg et al. 2005[ [Bisschop et al. 2006P show that N2 is 
only marginally more volatile than CO. This result throws 
doubt on much recent modelling of the chemistry of 
prestellar cores and poses the question of whether another 
more volatile and abundant form of nitrogen is responsi- 
ble for the survival of N-containing species to high density. 
Another possibility is that atomic nitrogen does not stick 
effectively to grain surfaces. In recent studies (Flower et al. 
2005, 2006), we have explored this suggestion and found 
that the best fit of the models to the observed abundances 
of NH3 and N2H+ is obtained when the sticking probabil- 
ities of both atomic N and O are low and, additionally, 
the mean grain surface area per H atom is a factor 5-10 
lower than in the diffuse interstellar medium; this implies 
that considerable grain growth has occurred prior to the 
prestellar-core phase. 
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In all of this discussion, the abundance of oxygen- 
containing species, such as OH and water, has been some- 
what neglected, and for good reason. OH and water are 
hard to observe with reasonable (arc min or better) reso- 
lution in cold clouds (see, for example. Goldsmith & Li 
2005, Bergin & Snell 2002). Thus we do not presently 
know whether oxygen-containing species survive in the 
gas phase in the regions where N2H+ and NH3 are ob- 
served. In the long run, observations of water with the 
HIFI instrument on the Herschel satellite or of OH with 
the Square Kilometer Array (SKA) may provide an an- 
swer. In the meantime, we need to consider alternative ob- 
servational strategies, using chemical models of prestcllar 
cores as a guide. If atomic oxygen sticks to grain surfaces 
no more efficiently than atomic nitrogen, OH and H2O 
might be expected to survive in the gas phase similarly to 
NH3 and N2H+. 

An important aspect of this problem involves NO, 
which is believed to be a key molecule in the nitrogen 
chemistry. NO mediates the transformation of atomic into 
molecular nitrogen, in the reactions 



N -I- OH ^ H + NO 
and 

N -I- NO ^ O + N2 



(1) 



(2) 



Atomic and molecular nitrogen are expected to be the ma- 
jor forms of elemental nitrogen in the gas phase. Under 
conditions in which reactions H]) and dominate the 
formation and destruction, respectively, of NO, the ratio 
of the abundance of NO to that of OH is given simply by 
the ratio of the rate coefficients for these reactions, fci/fc2; 
ki/k2 ~ 1, providing that there are no energy barriers to 
these reactions, which might reduce differentially the val- 
ues of their rate coefficients at low temperatures. Then, 
one might consider NO to be a proxy for OH, providing 
information on 0-containing species in the gas phase. OH 
forms through a sequence of reactions which is initiated 
by 0(H;^, H2)0H+ (and the analogous reactions with the 
deuterated isotopes, H2D+ and D2H+), followed by hydro- 
genation by H2 and dissociative recombination with elec- 
trons. OH is destroyed principally in the reactions 0(0H, 
02)H and N(OH, NO)H. Thus, the abundance of OH is 
tied to the abundances of Hj^, H2D+ and D2H+, and OH 
is a minor oxygen-containing species. 

NO is a known interstellar molecule, which has been 
detected in L183 with a fractional abundance of around 
lO"'^ (see Gerin et al. 1992, 1993). The value quoted by 
these authors was an overestimate, as it was based on 
the C^^O column density, and CO is now known to be de- 
pleted in L183 (Pagani et al. 2004, 2005). Nonetheless, the 
fractional abundance of NO is high, which encouraged us 
to investigate the spatial variation of the NO abundance 
in the cores LI 83 and LI 544. For these objects, estimates 
of the column densities of molecular hydrogen, A^(H2), 
which are independent of the degree of depletion, may be 
made from the existing maps of the dust emission. Our 



objective was to compare the spatial variations of the NO 
and H2 (from the dust) column densities. In addition, we 
observed simultaneously with NO the J = 1 — » tran- 
sitions of N2H+, which are known to follow roughly the 
dust emission. 

In Sect. 2, we describe our observations and, in Sect. 3, 
give the observational results. Comparisons with model 
predictions are made in Sect. 4, and our concluding re- 
marks are to be found in Sect. 5. 



2. Observations 

The observations were carried out during two sessions, in 
July and August 2005, using the IRAM 30 m telescope. 
We observed simultaneously with the facility 3 mm, 2 mm, 
and 1.3 mm receivers. The 3 mm receiver was tuned to the 
N2H+ J = 1 ^ 93.17632 GHz fine, the 2 mm receiver 
to the centre of the NO J = 3/2 1/2 H+ and H" 
transitions at 150.36146 GHz, and the 1.3 mm receiver to 
the centre of the NO J = 5/2 ^ 3/2 H+ and H" tran- 
sitions at 250.61664 GHz. The half-power beam widths 
at these three frequencies are 27 , 16 , and 10 , respec- 
tively, and the corresponding main-beam efficiencies are 
0.78, 0.68, and 0.48. Pointing was checked at intervals of 
roughly 2 hours and was found to be accurate within 2 . 
We used central positions (J2000) of R.A. = 05''04"16^9, 
Dec. = 25°10'47".7 for L1544 and R.A. = 15'>54'°08^8, 
Dec. = -02°52'44".0 for L183, and the offsets specified 
below are with respect to these positions. 

We used the facility SIS receivers in frequency- 
switching mode with system temperatures of 150-180 K at 
3 mm, 350-450 K at 2 mm, and 1000-2000 K at 1.3 mm. 
At 3 mm, we used the facility VESPA autocorrelator as 
a spectral backend with 20 kHz channel spacing over a 
20 MHz band, allowing us to cover all 7 hyperfine satel- 
lites of the N2H+(J =1^0) transition. At 2 mm and at 
1.3 mm, the NO lines are too widely separated to be cov- 
ered in one band with adequate resolution. We therefore 
split VESPA into 6 parts of 20 MHz bandwidth and 20 kHz 
resolution at 2 mm and 3 parts with 40 MHz bandwidth 
and 40 kHz resolution at 1.3 mm. The different sections of 
the autocorrelator were offset so as to cover the frequencies 
and transitions summarized in Table [T] The correspond- 
ing velocity resolutions were 0.04 km at 2 mm and 
0.05 km s-i at 1.3 mm. The data were processed using 
the GILD AS package^. 



3. Results 

In Appendix A, we describe the method used to deter- 
mine the column densities of N2H+, which are used in the 
chemical analysis presented in the following Section. Here, 
we discuss the derivation of the column densities of NO 
for each of the two prestellar cores which we observed. 



See URL http:/ /www.iram. fr /IRAMFR/GILD AS / 1 
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Table 1. NO line parameters 



Transition 

,f' ^j",f" 



Frequency 
GHz 



Line Strength^ 



H-* 



3 5 



2' 
3 



2 

_ 1 

2 ' 2 

H" 

3 5 13 

2 ' 2 2 ' 2 

3 111 
2 ' 2 2 ' 2 

n+ 

5 7 3 5 
2 ' 2 2 ' 2 

H" 

7 3 5 

2 2 ' 2 

5 3 3 

2 2 ' 2 



5 

2 ' 
5 

2 ' 



150.17646 
150.19876 
150.21874 
150.22565 

150.54646 
150.58055 

250.43684 

250.79643 
250.81561 



0.500 
0.186 
0.148 
0.148 

0.500 
0.148 

0.445 

0.445 
0.280 



(1) Gerin et aL (1992) 



3.1. L1544 

In the time available, rather than attempting to map the 
NO emission, we observed along two cuts, in the NW- 
SE and NE-SW directions, centred on the origin given in 
Sect, m these directions correspond roughly to the major 
and minor axes, respectively, of the dust emission, which 
we assume to be a good representation of the total col- 
umn density distribution. Furthermore, the intensities of 
the NO 150 GHz (2 mm) lines are expected to yield a good 
approximation to the true NO column density distribution 
(see IGerin et al. 1992p . as long as the lines are optically 
thin; this is likely to be the case, given their relative weak- 
ness (the low J levels are expected to be thermalized in 
L1544). Furthermore, the relative intensities of the dif- 
ferent hyperfinc components are observed to be close to 
the values expected on the basis of the line strengths in 
Table [TJ Averaging the corresponding transitions from the 
n+ and n~ bands, we find a ratio of 0.46 ± 0.07 between 
the J = 3/2 ^ 1/2 F = 3/2 ^ 1/2 and F = 5/2 ^ 3/2 
lines (at zero offset, (0,0): see Tabled]), as compared to 0.37 
using the line strengths in Table [1] Similarly, the ratio of 
the average of the three 2 mm lines with strength 0.148 in 
Table[T]to the two lines with strength 0.500 is measured to 
be 0.30 ±0.04, as compared with 0.30 in the optically thin 
LTE limit. Thus, wc find that the NO hyperfine satellite 
intensities are almost consistent, to within the error bars, 
with the low optical depth limit. Wc can then convert 
the integrated intensity I{2 mm) in the J = 3/2 ^1/2 
F = 5/2 ^ 3/2 (average of n+ and H") into an NO col- 
umn density for a given excitation temperature, rcx(NO). 

Because NO has a small dipole moment, the criti- 
cal density (at which the rates of coUisional and radia- 
tive de-excitation are equal) is of order lO** cm^'^, and 
it seems reasonable to suppose thermalization of the NO 
levels of interest to us here. Based on the NH3 results 

8.75 K and find 



Table 2. NO transitions observed towards L1544 (0,0), 
with uncertainties in parentheses. The linewidths of the 
150 GHz transitions were determined by means of a gaus- 
sian fit, using the CLASS procedure; the linewidths of 
the (noisier) 250 GHz lines were fixed to that of the 
150.17646 GHz transition. 



Frequency 


Integrated intensity 




Av 


GHz 


K km s'l 


km s~^ 


km s~^ 


150.17646 


0.22(0.02) 


7.18(0.01) 


0.38(0.02) 


150.19876 


0.09(0.01) 


7.14(0.02) 


0.35(0.04) 


150.21874 


0.05(0.01) 


7.12(0.02) 


0.30(0.04) 


150.22565 


0.06(0.01) 


7.12(0.02) 


0.26(0.04) 


150.54646 


0.20(0.01) 


7.11(0.01) 


0.34(0.02) 


150.58055 


0.04(0.01) 


7.13(0.03) 


0.34(0.06) 


250.43684 


0.26(0.04) 


7.28(0.05) 


0.38 


250.79643 


0.22(0.04) 


7.26(0.04) 


0.38 



of Tafalla et al. (2002) we assume T, 



that the column density of NO cm ^ in L1544 should be 
given approximately by 

iV(NO) = 4.3 X 10^^ I{2 mm) 

where I{2 mm) K km s^-^ denotes the intensity of the 
NO 2 mm (150 GHz) lines (see lGerin et al. 19921 table 3). 
We note that new VLA observations by Crapsi et al. 
(2005b) suggest that the temperature decreases towards 
the dust peak; but, as we shall see below, NO is likely to 
be depleted at high densities, and so we use the single- 
dish temperature estimate. An independent check of the 
NO excitation is provided by our measurements of the 
J = 5/2 ^ 3/2 250 GHz lines. From Tabled we conclude 
that /(1. 2 mm)//(2 mm) = 1.14±0.25, where, in addition 
to the formal errors given in the Table, we have assumed a 
15 percent calibration error. Here, 7(1.2 mm) is the aver- 
age of the integrated intensities of the two J = 5/2 3/2 
F = 7/2 5/2 transitions. Using the formulation of 
Gerin et al. (1992), we find that the corresponding ex- 
citation temperature between the J = 5/2 and J = 3/2 
levels is Tox = 13.5 ± 4 K, which is higher than the single- 
dish NH3 measurement but consistent with it, to within 
the combined error bars. 

In Fig. [U we show our observations of NO at 2 mm 
for both the NE-SW and NW-SE cuts, together with 
the corresponding measurements of the dust emission 
(from the map of [Ward-Thompson et al. 1999P and our 
3 mm observations of N2H+ J,Fi,F = 1,0,1 0,1,2 
(93.17632 GHz) transition. In essence, these data show 
that, while the dust and the N2H+ emission peak at the 
(0,0) position and have an angular size of about 1 by 2 
(0.04 pc by 0.08 pc at the adopted distance of 140 pc of 
L1544; [Crapsi et al. 2005aP , the NO emission extends over 
at least 3 and weakens, relative to the N2H"'" and dust 
emission, towards the (0,0) position. We can quantify this 
result by using estimates from Fig. [1] of the flux in the 
dust continuum emission. We derive a peak of emission at 
1.2 mm of 225 mJy in the beam, which corresponds, for 
an assumed dust absorption coefficient k per unit mass 
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Fig. 1. Plots of the average of the integrated intensities 
of the NO n+ and n J, F ^ 3/2,5/2 J,F = 1/2,3/2 
lines (see Table [ij, denoted by full curve joining circles 
with error bars, which have approximately the same size 
as the symbols, plotted against offset from the (0,0) posi- 
tion in L1544. Increasing offset is in the sense of increasing 
R.A. The top panel shows the cut along the NE-SW ma- 
jor axis, and the bottom panel the cut along the NW-SE 
minor axis. We show also for comparison the integrated in- 
tensity of the 93.17632 GHz transition of N2H+ (squares 
joined by long broken curve) and the 1.2 mm dust emis- 
sion (diamonds joined by short broken curve), from Ward- 
Thompson et al. (1999). The error bars plotted correspond 
to RMS deviations from the fits (gaussian for the lines); 
thus we neglect systematic errors arising from the calibra- 
tion, for example. The NO measurements involved long 
integration times, and the two components (11+ and H^) 
were averaged, whence the small error bars in this case. 
The intensities derive from main-beam brightness tem- 
peratures. 



Fig. 2. Plots of the average of the integrated intensities 
of the NO n+ and H" J, F = 3/2,5/2 ^ J, F = 1/2,3/2 
lines (see Table [T]), denoted by full curve joining circles 
with error bars, against offset from the (0,0) position in 
L183. Increasing offset is in the sense of increasing R.A. 
The cut is along the E-W (minor) axis. We show also for 
comparison the integrated intensity of the 93.17632 GHz 
transition of N2H+ (squares joined by long broken curve) 
and the 1.2 mm dust emission (diamonds joined by short 
broken curve), from Pagani et al. (2004). The intensities 
derive from main-beam brightness temperatures. 

density p of gas, n/p = 0.01 cm^ g~^, to an H2 column 
density of 1.2 x 10^^ cm~^; we have assumed that the 
dust temperature is equal to the excitation temperature, 
Tex = 8.75 K. This conversion factor was used to derive 
the column density ratio, a;(NO) = iV(NO)/7V(H2), shown 
in Fig. [3l We see that the observed column density ratio is 
approximately 8 x 10~^ at the dust peak and increases to 
values of 3 X 10~® to 5 x 10~^ at offsets of around 1 (equiv- 
alent to 0.04 pc, or 8000 AU) along the minor axis, where 
the H2 column density, derived from the dust emission, is 
a factor of 6 lower than at the peak. Along the major axis, 
the variation is less marked, but there is also an increase 
in the column density ratio away from the peak of the 
dust emission. We suggest in the next Section that this 
behaviour is probably related to the partial depletion, at 
high densities, of the main forms of oxygen. 

3.2. L183 

Our observed line parameters, towards the (0,0) position 
of L183, are given in Table [3] In this case, only one cut 
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Table 3. NO transitions observed towards L183 (0,0), 
with uncertainties in parentheses 



Frequency 


Integrated intensity 




At; 


GHz 


K km s"^ 


km s^^ 


km 


150.17646 


0.47(0.015) 


2.37(0.01) 


0.38(0.02) 


150.19876 


0.18(0.015) 


2.37(0.02) 


0.44(0.04) 


150.21874 


0.14(0.014) 


2.29(0.02) 


0.38(0.04) 


150.22565 


0.14(0.014) 


2.35(0.02) 


0.32(0.04) 


150.54646 


0.46(0.016) 


2.38(0.01) 


0.38(0.02) 


150.58055 


0.12(0.012) 


2.38(0.02) 


0.29(0.04) 


250.43684 


0.27(0.060) 


2.33(0.01) 


0.11(0.03) 


250.79643 


0.32(0.060) 


2.38(0.02) 


0.27(0.06) 



was possible in the time available, in an E-W direction, 
and the results are shown in Fig. [21 

We have used a procedure to derive the NO abundance 
similar to that for L1544. We assume a temperature of 9 K, 
based on the NH3 results of Ungerechts et al. (1980), in 
order to compute both the NO and H2 column densities. 
We use additionally the dust emission map of Pagani et 
al. (2004). As discussed below, our NO column density 
peak in L183 appears to be shifted by about 10 in R.A., 
relative to the dust peak (roughly 1000 AU for an assumed 
distance of 110 pc: [Franco 1989p . We infer column density 
ratios, iV(NO)/7V(H2), in the range 1 x 10"® to 3 x 10"^, 
similar to the values derived for L1544. 

4. Modelling the NO abundance distribution 

The simple models discussed below are based on the as- 
sumption of spherical symmetry, which is clearly a rough 
approximation for both of the prestellar cores which we 
consider. Indeed, Doty et al. (2005) have modelled L1544 
as a prolate spheroid, with an axis ratio of approximately 
2. Such departures from spherical symmetry imply that 
the timescale for collapse is different along each of the 
principal axes. Furthermore, determinations of the ther- 
mal pressure in these objects and estimates of their mag- 
netic field strengths suggest that the timescale for collapse 
is probably larger than for free-fall. Thus, the "free-fall" 
collapse model which we have adopted is undoubtedly a 
zero-order approximation. Nevertheless, for the purposes 
of comparing observations of different chemical species, 
and hence deducing the chemical and physical conditions 
in the cores, the present approach is probably adequate. In 
any case, to go further would necessitate a more complete 
knowledge of the dynamics of the cores. 

4.1. Empirical model 

First, we simulate the observed NO column density dis- 
tribution, relative to the II2 column density deduced from 
the dust emission (cf . Fig. [T]) , by means of a simple em- 
pirical model, according to which n(N0)/n(Il2) behaves 
as a step function, centred on the (0,0) position. The ra- 
tio a;(NO) = iV(N0)/iV(H2) is then computed, assum- 



ing spherical symmetry and convolving with the telescope 
beam, taken to have a width of 17 . In Fig. [3l we plot, for 
the case of L1544, two cases, one in which n(NO)/7i(H2) = 
5 X 10-9 for offsets r < 30" (4000 AU), and n(N0)/n(H2) 
= 3 X 10"* for r > 30 and the other in which the ra- 
tio n(N0)/n(H2) is assumed to be zero within the cen- 
tral zone. It may be seen from Fig. [3] that the value of 
n(N0)/n(Il2) = 5 X 10^9 for the abundance ratio in the 
central zone appears to be a reasonable upper limit, given 
the error bars on the observations. The asymmetries of 
the observed abundance ratio are not reproduced by the 
model, which assumes spherical symmetry. 

The central zone, in which a;(NO) decreases, is smaller 
than has been estimated for CO and some other species; 
the decrease occurs at a density n(H2) w 4 x 10^ cm"^. 
Assuming an exponential decrease of the fractional abun- 
dance with radial distance, Tafalla et al. (2002) derived 
characteristic (1/e) densities of 5.5 x 10^ cm"^ for C-'^^O 
and 1.7 X 10^ cm~^ for CS. Thus, NO appears to be in- 
termediate in its depletion characteristics between the C- 
containing species, which are strongly affected by CO de- 
pletion, and N-bearing species like N2H+ and NH3. 

An analogous model for our R.A. cut in LI 83, in which 
the Crapsi et al. (2005a) density distribution was adopted 
(a central density of 10^ cm~'^ and a power-law decrease 
in the density beyond a radius of 4200 AU), shows the 
displacement, to which we alluded above, of the observed 
minimum of the NO abundance ratio from the dust peak. 
In this case also, it will be necessary to map fully the NO 
column density distribution to take the analysis further. 

4.2. Chemical considerations 

In prestellar cores, the nitrogen chemistry depends on 
gas-phase conversion of atomic into molecular nitrogen. 
As was mentioned in the Introduction, the most efhcient 
scheme for this conversion is believed to involve NO as an 
intermediate, through reactions ([T]) and ^ above. The ra- 
tio of the abundance of NO to that of OH is then given by 
ki/k2, the ratio of the rate coefficients for these reactions. 

This relation no longer applies when the atomic ni- 
trogen abundance falls below the value for which the de- 
struction of NO by N, in reaction has the same rate 
as direct depletion of NO on to the grains, i.e. when 

nHfcdcp(NO) = n(N)fc2, 

or 

x(N) EE n(N)/nH = fcdcp(N0)/A;2. 
In these expressions, 

A:dcp(NO) = (ng/nH)7ragWth, 

where rig is the grain number density, Og is the grain ra- 
dius, and ?;th is the thermal speed of the NO molecules. 
Taking T== 10 K, we find A:dep(NO) 1.4 x 10-i*(0.5/ag), 
where Og is in /im and /cdep in cm"^ s^^. With k2 — 
3.5 x 10-" cm^ s-i (Baulch et al. 2005), we find that, 
when the fractional abundance of atomic nitrogen exceeds 
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10- 



o 
z 



lo-" 




(Sims 2005), and so the values of ki and k2 remain un- 
certain. As will be seen in Sect. 14.31 our observations of 
NO and N2H+ support the existence of a small barrier to 
reaction ^ , which inhibits the destruction of NO by N at 
temperatures as low as T = 10 K. 

We recall that our observations indicate that, unlike 
NH3 and N2H+, NO does not show a constant or even 
an increasing abundance, relative to H2, for densities 
3 X 10^ cm~'^. Instead, NO appears to decrease in 
relative abundance approaching the dust emission peak. 
Here we review some possible chemical explanations of 
this discrepancy. 



-100 -50 50 

Offset (Arc sec) 



100 




50 

Offset (Arc sec) 
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Fig. 3. Observed and predicted NO to H2 column den- 
sity ratio, a;(NO) A^(N0)/iV(H2), for a model of 
L1544 in which the Tafalla et al. (2002) density dis- 
tribution (cf. Sect. 14. 3p is adopted. Two different as- 
sumptions (indicated on the Figure) are made regard- 
ing the value of n(N0)/n(H2) in the central region, off- 
sets r < 30" (4000 AU): n(N0)/n(H2) = 5 IQ-^ (broken 
curve); n(N0)/n(H2) = (full curve). Outside this region, 
n(N0)/n(H2) = 3 10~* is assumed. The filled circles with 
error bars denote the observed values. 



a;(N) = 4 X 10"^, gas-phase destruction of NO by N domi- 
nates, whereas, below this value, NO is removed predomi- 
nantly by freeze-out on to the grains. A similar argument 
can be applied to OH, and so the abundances of NO and 
OH remain related in both regimes. 

The rate coefficients for both of the reactions ([T|) and 
(HJ must be sufficiently large at low temperatures to ensure 
the conversion of N into N2, the precursor of N2H+ and 
NH3; but the measurement of rate coefficients for neutral- 
neutral reactions at very low temperatures is demanding 



OH may be under-abundant in the high density 
regions of at least some prestellar cores. This sit- 
uation could arise for the same reasons that O2 
and water are much less abundant than expected in 
prestellar cores (see, for example, [Pagani et al. 20031 
ISnell et al. 2000p . Surplus oxygen - that is, oxygen 
which is not bound in CO - may have frozen out prior 
to the start of the collapse, giving rise to a situation 
whereby, once CO has also frozen out, there is more 
N than O in the gas phase; this assumes that N does 
not freeze out simultaneously with CO. In the lower 
density regions, in which CO is still present in the gas 
phase, the C/0 elemental abundance ratio would be 
close to unity, and the chemistry would have "carbon- 
rich" characteristics, with, for example, large abun- 
dances of cyanopolyynes and similar species. In this 
case, the main source of atomic oxygen is the reaction 
of CO with Hc+ , which is produced by cosmic ray ion- 
ization of He, and hence the fractional abundances of 
both OH and NO are dependent on the cosmic ray flux. 
Also, their abundances are determined to some extent 
by how rapidly atomic oxygen sticks to the grains, an 
issue which will be discussed in Sect. 14.31 below. 
There may be chemical pathways from N to N2 which 
we (and others) have overlooked. Any conversion of N 
into N^ , through charge exchange with He^ , for exam- 
ple, would lead to the production of NH and other hy- 
drides and the subsequent formation of N2 in neutral- 
neutral reactions. However, trial calculations which we 
have carried out suggest that this particular route is 
not likely to be significant. More promising is the pos- 
sibility that, under conditions in which the gas-phase 
abundance of elemental carbon approaches that of el- 
emental oxygen, the reactions 



N + CH ^ H + CN 
and 

CN ^ C 



N 



No 



(3) 



(4) 



assume importance relative to reactions ([T]) and 
We consider this possibility further in Sect. 14.31 



We conclude that the main forms of oxygen (like the 
main forms of carbon) appear to have frozen on to dust- 
grain surfaces at densities approaching 10^ cm^'^ and. 
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moreover, that some form of nitrogen is more volatile and 
hence more abundant in the gas phase than any form of 
oxygen. Another conclusion is that, in the densest parts of 
prestellar cores, oxygen exists mainly as solid H2O, CO, 
and perhaps CO2, with a fractional abundance in the solid 
state which is larger than the values derived for lower den- 
sity 3 X 10^ cm~^) material. Having eliminated other 
possibilities, we arrive at the proposal that atomic nitro- 
gen does not stick effectively to icy dust grains, even at 
kinetic temperatures close to 10 K (although we note that 
a small abundance of OH or CH is also necessary to the 
formation of N2). This proposal is not in conflict with 
the recent laboratory studies of Bisschop et al. (2006) and 
Oberg et al. (2005), which suggest that N2 freezes out 
similarly to CO. Then, atomic N is likely to be the main 
repository of nitrogen in the gas phase. 



4.3. Free-fall collapse model 

To test some of the above ideas, we have carried out 
calculations of the variation of NO, N2H+, and NH3 with 
density in the framework of the free-fall collapse model 
presented by Flower et al. (2006). In this study, it was 
found that a good fit to the observations of LI 544 in NH3, 
N2H+, and their deuterated isotopes required a grain sur- 
face area per H-atom roughly an order of magnitude less 
than in the diffuse interstellar medium. These simulations 
assumed a single grain size (radius) of typically 0.5 fim. 
It was assumed also that the sticking coefficients for both 
atomic nitrogen and oxygen are significantly less than 
unity. Then, NH3 and N2H+ remain in the gas phase at 
densities for which CO and other C-containing molecules 
are already depleted. In the present Section, we use the 
model developed by Flower et al. (2006) and examine its 
predictions for the NO abundance. When so doing, we 
take account of the fact that our observations indicate 
that NO is intermediate between the C-containing 
species and molecules like NH3. The complete chem- 
istry file and list of chemical species are available from 
|http : / / massey.dur.ac.uk/drf/ protostellar/species-chemistry. 

One possibility which has been investigated is that our 
assumed initial abundances (calculated in steady state for 
a density nn — 10'* cm~'^) influence the results signif- 
icantly. In particular, we have compared our 'standard 
model', in which the ratio of elemental carbon to oxygen 
in the gas-phase is initially 0.67, to a model in which this 
ratio approaches 1. In this latter case, the underlying as- 
sumption is that oxygen which is not initially in the form 
of gas-phase CO is in water ice. The results of these cal- 
culations are shown in Fig. [H We see from this Figure 
that the variation of the fractional abundance of NO, 
n(NO)/riH, with nn is very similar in the two cases; this is 
because NO forms in the reaction ([T]) of N with OH. The 
OH radical is formed and destroyed in binary reactions in- 
volving atomic O, and so its fractional abundance is insen- 
sitive to the total amount of oxygen in the gas phase. On 
the other hand, there is a large decrease in the fractional 



abundance of O2 when the ratio nc/no is increased. The 
computed fractional abundance of O2 becomes just about 
compatible with the upper limit of n(02)/n(H2) 10^^ 
derived by Pagani et al. (2003), from observations of a 
number of prestellar cores, in the 'carbon-rich' case. For 
this reason, we adopted nc/no — 0.97 « 1 initially, in the 
gas phase, in subsequent models. 

We have already mentioned in Sect. 14.21 that N2 is 
the precursor of N2H+. Therefore, the initial value of 
the fractional abundance of N2H+ depends on the ini- 
tial fractional abundance of N2. Because reactions ([T])- 
([4]), which transform N into N2, involve minor species, 
OH and CH, this transformation process is slow, requir- 
ing approximately 5 x 10^ yr to attain equilibrium at a gas 
density nn = 10^ cm~'^. Thus, it is far from clear that the 
ratio n(N)/n(N2) can be assumed to be initially in equi- 
librium (cf. Flower et al. 2006). Accordingly, we have var- 
ied the initial value of this ratio from n(N)/n(N2) << 1, 
obtained assuming that the chemistry has reached equi- 
librium, to n(N)/n(N2) >> 1, with a view to improving 
the agreement with the fractional abundance of N2H+ de- 
duced from our observations. 

The observed quantities are the column densities of 
NO, N2H+ and H2, as functions of position. In order to 
compare the free-fall model with the observations, we re- 
late the offset position to the gas density by means of the 
relation 

n(0) 
= YT(^ 

(Tafalla et al. 2002) , where n is the gas density, r is the off- 
set from centre, r — Q, n(0) is the central density of molec- 
ular hydrogen, and r2 is the radial distance over which the 
density decreases by a factor of 2, relative to the central 
density. From the 1.2 mm dust continuum emission map of 
L1544, Tafalla et al. derived r2 = 20" (equivalent to 0.014 
pc at the distance of L1544), a = 2.5, and n(0) = 1.4 x 10^ 
cm~'^. In view of the freeze-out density indicated by the 
free-fall model (Fig. [5]), we have adopted a central density 
T?.(Ol — 0.5 X 10^ cm~^ of molecular hydrogen, equivalent 
-fi^-A&'(0) = 1.0 X 10^ cm-3. We consider that the cen- 
tral density is uncertain to at least this factor of 2.8 (cf. 
Bacmann et al. 2000, Crapsi et al. 2005a, Doty et al. 2005). 
We note also that the density n{r) which derives from the 
above expression is the maximum value along the line of 
sight at a given offset, r, from centre. 

From our earlier discussion (Sect. 14. 2p . we see that 
NO might be sensitive also to the sticking coefficient 
of atomic oxygen, 5(0), and so we compare results for 
5(0) = 0.1 and S{0) = 1.0 in Fig. [3 The results of 
the calculations are, indeed, sensitive to S{0): NO van- 
ishes from the gas phase at a density of roughly 10^ cm^'^ 
when ^(O) = 1.0 but survives to at least 3 times that 
density when S{0) — 0.1. Of course, CO is unaffected 
by the oxygen sticking probability. In these calculations, 
S{C) = S'(N) = 0.1 and S = 1.0 for all other neutral 
species. The model is 'carbon-rich' {nc/no = 0.97 ini- 
tially in the gas phase), and hence the route to N2, through 
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C/O = 0.67 




c/o = i 




Fig. 4. Results of models of the chemistry in a core collapsing on the free-fall timescale. Results are shown for 
models in which (a) nc/no — 0.67 and (b) nc/no — 0.97 « 1 initially. In both cases, the sticking probabilities 
^(N) = S{0) = 0.1 and S* = 1.0 for all other neutral species. 



the reactions of CH and CN with N (reactions ^ and 
(0])), assumes importance, relative to the reactions of OH 
and NO with N (reactions ^ and The initial abun- 
dance of CH exceeds that of OH by an order of magni- 
tude. It is noticeable that CH and CN are considerably 
more abundant at densities of the order of 10^ cm""^ in 
the models in which 5(0) = 1.0; atomic oxygen destroys 
these species, forming CO. A consequence is that, if S{0) 
is large, CN may be expected to be depleted less than 
other C-containing species towards the density peaks of 
prestellar cores. Finally, we see that the close link between 
the fractional abundances of NO and OH, mentioned in 
the Introduction and discussed in Sect. 14. 2[ is apparent in 
Fig. El 

In Fig. [S] are compared the observed values of the col- 
umn density ratios A^(NO)/A'^h and iV(N2H+)/7VH, where 
Nyi — 2A^(H2), with the computed values of the fractional 



abundances n(NO)/riH and n(N2H+)/nH. In this Figure, 
S{C) = S'(N) = 5(0) = 0.1. The horizontal error bars rep- 
resent a factor of 3 uncertainty in nn, whilst the vertical 
error bars incorporate the observational uncertainties in 
the line intensity or flux measurements, together with an 
assumed 15 percent calibration error. The observed points 
fall on two branches, corresponding to opposite sides of 
the central position. The emission from L1544 is observed 
to be asymmetric about the centre, whereas the empirical 
model of the density distribution, given above, assumes 
radial symmetry. (The lower branch corresponds to neg- 
ative offsets (in R.A.) for both NO and N2H+.) It may 
be seen that increasing the initial value of the atomic to 
molecular nitrogen ratio has the effect of yielding much 
improved agreement between the observed and computed 
profiles of N2H+. On the other hand, the profile of NO 
is displaced downwards, owing to the enhanced rate of its 
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S(C) = S(N) = S(0) = 0.1 S(C) = S(N) = S(0) = 0.1 





Fig. 5. Results of models of the chemistry in a core collapsing on the free-fall timescale. Results are shown for models 
in which the grain sticking coefficient of atomic oxygen is set equal to S{0) — 0.1 (upper panels) and S{0) = 1.0 
(lower panels). In both cases, S{C) — S{N) = 0.1 and S — 1.0 for all other neutral species; nc/no — 0.97 « 1 initially 
in the gas phase. 



destruction by N, and the computed profile now falls well 
below the observations. However, if a small barrier intro- 
duced into the rate coefScient for the reaction of NO with 
N, the computed NO profile becomes once again compati- 
ble with the observations, as may be seen in Fig. [6l In the 
calculations presented in this and the following Figures, we 
have adopted the rate coefficient for the reaction of NO 
with N which is given in the current version of the UMIST 
chemical reaction database (LeTeuff et al. 2000), rate05, 
namely ^2 = 3.75 x 10""exp(-26/r) cm^ s~\ which de- 
rives from the study of Duff & Sharma (1996). This small 
reaction barrier is sufficient to inhibit the destruction of 
NO by N significantly at T = 10 K. We note that it is 
presently very difficult to confirm or exclude the existence 
of such a small barrier by either laboratory measurements, 
which extend to insufficiently low temperatures, or theo- 



retical calculations, which would need to determine the 
transition state energy to an accuracy of 100 /ih. 

In Fig. [7l we investigate the effect of varying the grain 
sticking coefficient of atomic oxygen, in the range in the 
range 0.1 < S{0) < 1.0. Once again, S(C) = S{N) = 0.1. 
We see that the general trends of the observations are 
for 7V(NO)/7Vh to increase with decreasing density (in- 
creasing distance from the centre), whereas A^(N2H"'")/A^h 
remains approximately constant or, possibly, tends to de- 
crease with decreasing density; this is in accord with the 
comments made in the previous Section. The models in 
which S{0) > 0.1 predict that the fractional abundance 
of N2H+ is roughly constant, which, given the large er- 
ror bars, is in reasonable accord with the observations. As 
S{0) increases towards unity, the amount of oxygen re- 
maining in the gas phase at high densities decreases and 



10 



M. Akyilmaz et al.: Depletion of NO 





Fig. 6. A comparison of the column densities of NO and N2H+, relative to Nh ~ 2N{H2), derived from observations 
along the minor axis of L1544 (points with error bars), with the fractional abundance of NO, n(NO)/nH, and of 
N2H+, n(N2H+)/nH, predicted by free-fall models in which the initial value of the ratio n(N)/n(N2) is varied from 
its equilibrium value, n(N)/n(N2) = 1/60 (dashed curve) to n(N)/n(N2) = 18 (long-dashed curve). The effect of 
introducing a small (26 K) barrier in the rate coefficient of the reaction of NO with N (reaction ^) is also illustrated 
(full curve); the introduction of the barrier has no perceptible effect on N2H+, where the long-dashed and full curves 
overlap. In these calculations, 5'(C) = 'S'(N) = S{0) = 0.1 and S = 1.0 for all other neutral species; initially, 
n-c/iT'O = 0-97 « 1 in the gas phase. Note that the abscissa has been reversed to reflect the decrease of the density 
with increasing distance from the centre. 
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the rate of formation of NO, in reactions H]) and de- 
creases also. When S{0) = 0.5, n(NO)/nH = 3.5 x IQ-^ 
when nil = 10^ cm"^, and n(NO)/nH = 2.9 x 10^^ when 
riH = 10^ cm^'^; these values are consistent with those 
deduced from the empirical model (Sect. 14. ip . However, 
we wish to emphasize that we do not consider ourselves 
in a position to "determine" the sticking coefficients: this 
should be done in the laboratory. We have treated the 
sticking coefficients of C, N, and O simply as parameters 
of the model, varying their values in an attempt at fitting 
the available observations. 

In our recent study (Flower et al. 2006), we found 
that satisfactory agreement could be obtained with the 
observed high levels of deuteration of the nitrogen- 
containing species NH3 and N2H+, providing that the 
grains were sufficiently large (cg « 0.50 ^m, as in the 
present calculations) and that the sticking probabilities 
of atomic nitrogen and atomic oxygen were both less than 
unity {S <^ 0.3). For the reasons given above, we have since 
modified the initial values of the C/0 and N/N2 ratios, as 
well as reducing the sticking probability of atomic carbon 
to 5(C) = 0.1. We have seen (Fig. that the observed 
variations of the column densities of NO and N2H+ can 
then be reproduced, but it is necessary to estabhsh that 
this success is not at the expense of the earlier agreement 
with the degrees of deuteration of NH3 and N2H+. 

In Fig. [8] are plotted the relative abundances of the 
deuterated and non-deuterated forms of NH3 and N2H+ 
which derive from the present calculations, adopting 
nc/no = 0.97 and n(N)/n(N2) = 18 initially, together 
with 5(C) = 5(N) = 0.1; the sticking probability of 
atomic oxygen is varied in the range 0.1 < 5(0) < 0.5. 
The observational points and their error bars are identical 
to those adopted in our earher paper (Flower et a. 2006, 
figure 5c). It may be seen from Fig. [8] that satisfactory 
agreement with the observed levels of deuteration is still 
obtained when 5(0) « 0.3. This Figure emphasizes the 
fact that the degree of deuteration of NII3 increases with 
increasing gas density. The fractional abundance of ND3 , 
for example, has a distinct maximum at densities of order 
10^ cm~^. As kinematical parameters, such as the collapse 
speed, might be expected to vary with radial distance (and 
hence density) in the core, the widths of emission lines 
from successive stages of deuteration of NH3 could pro- 
vide information on the variation of the infall speed with 
the gas density. 

We note that the effect of subsequent coagulation on 
the size of initially large grains - specifically, Og = 0.50 /im 
in the present models - is small, over the range of com- 
pression considered (10^ < 7^H < 10^ cm~^), and has little 
consequence for the freeze-out timescale. 

5. Concluding remarks 

Our principal observational conclusion is that NO, unlike 
NH3 and N2II+, is depleted towards the density peaks 
of the prestellar cores L1544 and L183. This result was 
surprising, in that NO was believed to be an important 



intermediate of the nitrogen chemistry. In particular, the 
path from atomic to molecular nitrogen was believed to 
involve both OH and NO. Our results cast some doubt 
on this assumption, suggesting that molecular nitrogen 
may form from CH and CN at high densities. If so, CN 
is an important intermediate in the nitrogen chemistry. 
Further mapping of CN in prestellar cores, extending the 
observations of L183 by Dickens et al. (2000), would be 
extremely useful. Our models suggest a fractional abun- 
dance n(CN)/nH ~ 10"^*^ towards the peak of L1544, cor- 
responding to a column density of the order of 10^^ cm~^. 
In the case of L183, the central depression in the fractional 
abundance of NO is offset from the dust peak. This object 
needs further study: complete maps of the NO distribution 
are required. 

In an attempt to explain our observations within the 
framework of a free-fall collapse model, we have varied 
the values of the sticking probabilities of atomic C, N and 
O, and also the initial values of the elemental C/0 and 
the atomic to molecular nitrogen abundance ratios in the 
gas phase. We find that 

— the prestellar core L1544 is likely to be 'carbon-rich', 
in the sense that nc /f^o ~ 1 in the gas phase. An initial 
nc/no = 0.97 leads to values of n(02)/n(H2) which 
are in much better agreement with the observed upper 
limits than when the gas is oxygen-rich {nc/no = 
0.67). 

In objects in which the gas-phase C/0 abundance ra- 
tio is close to unity, the fractional abundances of C- 
rich species, such as the cyanopolynes, are expected 
to be enhanced relative to oxygen-containing species, 
such as O2 and H2O. Of the two cores which we 
have studied, this is true of L1544 (cf. Walmsley et 
al. 1980), whereas L183 is relatively carbon-poor (cf. 
Swade 1989). 

— The ratio n(N)/n(N2) at the start of collapse is much 
more than its value in equilibrium. This result is per- 
haps not surprising, given the long timescale (« 5 x 
10^ yr) required for the nitrogen chemistry to reach 
equilibrium at densities nn « lO'* cm~^. We note that 
Maret et al. (2006) have recently and independently ar- 
rived at essentially the same conclusion, viz. that most 
of the nitrogen in the gas phase is atomic, also from 
an analysis of observations of N2H+, in their case in 
the prestellar core B68 (but with implications for the 
deficiency of molecular nitrogen in cometary ices). If 
n(N)/n(N2) >> 1 in the gas phase, then our obser- 
vations of NO imply the existence of a small barrier 
(« 25 K) to the reaction of NO with N (equation 
Experimental demonstration of the existence of such 
a small barrier to this and, possibly, the other key 
neutral-neutral reactions H]), ([3]) and (|4]), whilst diffi- 
cult to achieve, is highly desirable. 

— The grain sticking probabilities of atomic C, N and, 
probably, O are significantly smaller than unity; val- 
ues of 5(C) = 5(N) = 0.1 and 0.1 < 5(0) < 1.0 were 
adopted in the calculations reported above. However, 
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the assumed values of the sticking probabilities for the 
atomic species are ad hoc and require confirmation by 
means of experimental measurements. Such measure- 
ments should presumably be made on CO ice, as in 
the experiments of Bisschop et al. (2006). 

Appendix A: Determination of the column 
densities 

The method for deriving the column density of NO, from 
observations of the hyperfine components of rotational 
transitions, assumed optically thin, has been presented by 
Gerin et al. (1992); we followed their method in the present 
paper. The case of N2H+ was considered by Caselli et 
al. (2002), who gave expressions for the total column den- 
sity of N2H+, for the case of a finite optical depth in the 
line and also in the optically thin limit, in their Appendix 
A. If the optical depth in a line of central frequency v and 
wavelength A is finite, then the total column density of the 
molecular ion is related to to the optical depth r at the 
line centre by 



agree to within about 30 percent with the results of more 
sophisticated analyses of the emission in the rotational 
multiplet, which allow for optical depth effects; this is the 
case even at the peak of the emission, where the assump- 
tion that the 1, 1, —> 0, 1, 1 transition is optically thin is 
least valid. Crapsi et al. (2005a) estimated the "total" op- 
tical depth (i.e. summed over all seven components) of the 
N2H+ J = 1 ^ multiplet to be r = 12.6 ± 0.7 towards 
the peak of L1544, which corresponds to an optical depth 
in the 1, 1, ^ 0, 1, 1 transition of 12.6/27 = 0.47 and an 
escape probability of 0.80 (and hence a 25 percent under- 
estimation of the column density owing to the "optically 
thin" assumption). They estimated, also from their fit, an 
excitation temperature Tex = 5 K for the J = 1 — s- 
multiplet in N2H"'". 

The hyperfine splitting of the rotational energy levels 
of ^''N2H+ arises from the interaction of the spin angular 
momentum, I, of the nitrogen nuclei and the rotational 
angular momentum, J, of the molecule. The angular mo- 
mentum coupling scheme is 
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where Aw is the full width of the line at half-maximum 
intensity (the line profile is assumed to be gaussian) , is 
the degeneracy (statistical weight) of the upper level of the 
observed transition, A(u — > 1) is the spontaneous radiative 
transition probability to the lower level, whose excitation 
energy is £'1; Q is the partition function, considered be- 
low. We have assumed implicitly a Boltzmann distribution 
of population at an excitation temperature (to be deter- 
mined) Tqx. In the limit of low optical depth at the line 
centre, the expression for the total column density is 



Fi 



N 



47rAi 



f— V 

Vln2/ 



Q 



mb 



1 



5u^(u 1) [eyiY>{hv / kB,Tc: 



1 



1 - exp(-/iz^/fcBTox) exp(-i?i/fcBTcx) 



(A.2) 



where Tmb is the main beam brightness temperature and 
Tbg is the temperature of the background radiation field, 
usually taken to be a black-body at Tbg — 2.73 K. 

The line width. At), and the brightness temperature, 
?nib, refer to a given component of the J = 1 — s- multi- 
plet, which is observed towards L1544 and L183 to com- 
prise seven resolved components. The full set of transitions 
is listed in Table lA.lj along with the corresponding line 
strengths. A comparison of the relative intensities of the 
components with the values expected in the optically thin 
limit (when they are given by the line strengths) shows 
that some of the lines are optically thick towards the dust 
emission peaks of these two sources. Accordingly, we have 
used the weakest transition, J,Fi,F = 1,1,0 0,1,1 of 
the multiplet, which should also be the most optically thin, 
to determine the column density of N2H"'". Previous work 
has shown that N2H+ column densities derived in this way 



The nuclear spin angular momentum quantum number is 
Ii = I2 — 1- When J > 1, there are three possible val- 
ues of Fi, namely Fi = J — 1, J, J + 1. On the other 
hand, when J ~ Q, Fi = Ii = 1. The possible values of 
F, the total angular momentum, are determined by the 
(weaker) coupling to the spin of the second (inner) nitro- 
gen nucleus. The quantum numbers J,Fi,F which iden- 
tify the hyperfine levels of the J = 1 and J — rotational 
states are listed in Table ETTl The optically allowed tran- 
' sirions are determined by the electric dipole selection rules 
e5^(^iz/%2jfej)(^JJ^^ = 0, when T = 0, is not allowed) 
governing the change in F in the transition. The hyperfine 
splitting of J = 0, which arises from the nuclear spin-spin 
interaction, is very small: see Green et al. (1974), Caselli et 
al. (1995). Consequently, allowed transitions from a given 
upper level, J — l,Fi,_F, are unresolved, and there re- 
main seven resolved hyperfine transitions between J = 1 
and J — 0, one from each of the hyperfine levels of J = 1. 

The partition function, Q, which appears in equa- 
tion (|A.2|) . is defined by 



gi exp 



E, 



(A.3) 



where the summation extends over all energy levels, and 
gi = 2Fi + 1. As Ei is approximately the same for all 
hyperfine levels belonging to the same J, the summation 
over gi for each value of J may be carried out, yielding 



Q = 9 ^ 5j exp 



Ej 



(A.4) 
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Table A.l. The J' — 1 J — hyperfine transitions 
at 93.2 GHz. The seven resolved lines are ordered in in- 
creasing frequency, from bottom to top. For each tran- 
sition, the quantum numbers identifying the upper level 
(with primes) and lower level are given, along with the 
corresponding line strength. The line strengths have been 
summed over the allowed transitions from a given hyper- 
fine level of J = 1, given that the hyperfine levels of J = 
are almost degenerate. 



Transition 
J'FiF' JFiF 



Line strength 
S{FlF',FiF) 



101 012 

121 Oil 
123 012 

122 Oil 

111 -> 010 

112 012 
110 Oil 



1/9 
1/9 
7/27 
5/27 
1/9 
5/27 
1/27 



where gj = 2J+l. The factor of 9 appearing in (jA.4p is 
the product of the degeneracies, 2/ + 1 = 3, arising from 
the spin of the nitrogen nuclei. It may be verified that, for 
a given value of J, 



^(2F, -1-1) ==9(2J-H) 



(A.5) 



Assuming that the hyperfine splitting is small com- 
pared with the rotational splitting of the energy levels 
(which is verified, in practice), the relationship between 
the probability of a given hyperfine transition within a ro- 
tational multiplet and the total probability of a transition 
between the rotational levels J', J is 

gF'A{ J'F[F' JFiF) = 9gj>A{ J' ^ J)S{F[F' , FiF){A.6l 

where gp' = 2i^' -I- 1 and gj' — 2 J' + 1; 9gj' is the to- 
tal degeneracy of rotational state J' (cf. equation (jA.SP ). 
S{F{F' , FiF) is the line strength, which determines the 
relative intensities of the hyperfine components of the ro- 
tational multiplet. The sum of the line strengths, over all 
components of the multiplet, is equal to unity. Because the 
hyperfine levels of J = are almost degenerate, the line 
strengths may be summed over over all allowed transitions 
from a given hyperfine level of J = 1. In Table ETTl the re- 
sulting line strengths are listed for the seven resolved tran- 
sitions from J' = 1 to J = in N2H+; they may be seen to 
be given by the ratio (2i^'-|-l)/ ^(2F'-|-1), where the sum- 
mation extends over the hyperfine levels of the upper rota- 
tional state, J'. When J' = 1, ^J2F^ + 1) = 27, as may be 



seen from equation (jA.SP : see Townes & Shawlow (1955) 
and Tatum (1986) for more detailed discussions of cou- 
pling schemes and line strengths. The specific case of 
N2H+ is considered in a recent publication of Daniel et 
al. (2006) (but note that J' in their equation (1) for the 
transition probability should read J). 



The derivation of the total column density of a 
molecule from observations of individual transitions - hy- 
perfine transitions between rotational states in the cases 
of NO and N2H+ - relies on correcting for the popula- 
tions of unobserved levels. When this is done via equa- 
tion (jA.3|) for the partition function, Q, it is assumed 
that the relative level populations may be characterized 
by a Boltzmann distribution at the excitation tempera- 
ture, Tox; this is the case when collisional de-excitation of 
the upper to the lower rotational state dominates spon- 
taneous radiative decay. The 'critical density' is defined 
as the density of the collisional perturber (predominantly 
H2 in the case of rotational transitions) for which the ra- 
diative and collisional decay rates become equal (once a 
correction has been made for the photon escape probabil- 
ity). For perturber densities higher than the critical den- 
sity, a Boltzmann distribution of the relative level popu- 
lations is approached, at the kinetic temperature of the 
gas. Corrections can be made for optical depth effects by 
means of either "large velocity gradient" programs (e.g. 
Crapsi et al. 2004) or more sophisticated (Monte-Carlo) 
radiative transport models (e.g. Tafalla et al. 2002). 

The peak column density of N2H+ in L1544 derived 
by Crapsi et al. (2005a), A^(N2H+) = (1.83 ± 0.19) x 
10^^ cm~^, who used the weakest component and the op- 
tically thin approximation, agrees to within the error bar 
with the determination of CaseUi et al. (2002), iV(N2H+) 
= 2.0 X 10^'^ cm~^, who used all seven hyperfine compo- 
nents and corrected for the finite optical depths. Thus, 
a simple approach to the determination of the column 
density of N2H+ can yield satisfactory results. To our 
knowledge, all line transfer models to date have assumed 
that the density structure of L1544 can be represented by 
a one-dimensional, spherically symmetric distribution. It 
would be interesting to compare with the predictions of 
radiative line transfer models which go beyond the ap- 
proximation of a single dimension, even though the angu- 
lar resolution and completeness of the currently available 
maps of the molecular emission would constrain the com- 
parison of such models with the observations. 

We are not aware of calculations of rate coefficients for 
rotational transitions in NO, induced by II2. However, cal- 
culations have been performed for OH - H2 (Offer & van 
Dishoeck 1992), which shares many physical characteris- 
tics with NO - H2, not least in the rotational structure of 
the target molecule. In the case of collisions with para-H2, 
which is the main form of H2 at low temperatures, the cal- 
culations of Offer & van Dishoeck (1992) indicate that the 
rate coefficient for rotational de-excitation at T = 50 K is 
of the order of 10^^^ cm'^ s~^ and not strongly dependent 
on the kinetic temperature, T. The rate of spontaneous ra- 
diative decay in the transitions of NO which we observed 
(cf. Table[I]) is of the order of lO"'^ s'^ (Gerin et al. 1992), 
which implies a critical density of H2 of the order of 10** 
cm~'^. In the regions of the protostellar cores considered 
above, n(H2) ^ 10^ cm~^, extending up to values in ex- 
cess of 10^ cm"'^. We conclude that the assumption of 
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Boltzmann equilibrium at the local kinetic temperature is 
justified in the case of NO. 

Consider now N2H+. Daniel et al. (2005) have recently 
computed rate coefficients for the rotational de-excitation 
of N2H+ by He; in the absence of data specifically for H2, 
we adopt their results as a guide. Their calculations indi- 
cate a value of 10~^^ cm^ s~^ for rotational transitions in 
N2H+, induced by He, at low temperatures, similar to the 
value above for NO - H2. On the other hand, the rate of 
spontaneous decay of J = 1 to J = is of the order of 
10"'"' (Dickinson & Flower 1981), 2 orders of magni- 
tude larger than for NO (N2H+ has a much larger dipole 
moment). Consequently, the critical density for N2H+ is 
of the order of 10® cm^'"*, although line trapping will tend 
to reduce this value towards the peaks of L1544 and L183. 
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Fig. 7. A comparison of the column densities of NO and N2H+, relative to Nh ~ 2N{H2), derived from observations 
along the minor axis of L1544 (points with error bars), with the fractional abundance of NO, n(NO)/nH, and of N2H+, 
n(N2H"'")/nH, predicted by free-fall models in which the grain sticking coefficient of atomic oxygen is varied: S{0) = 
0.1, 0.3, 0.5, 1.0. In all cases, S{G) — S'(N) ~ 0.1 and 5* = 1.0 for all other neutral species; initially, nc/no — 0.97 « 1 
in the gas phase, and n(N)/n(N2) = 18. Note that the abscissa has been reversed to reflect the decrease of the density 
with increasing distance from the centre. 
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Fig. 8. The abundance ratios of the deuterated forms of N2H+ and NH3, as functions of rin, together with the values 
observed in the prestellar core L1544 (cf. Flower et al. 2006, figure 5c). The grain sticking probability of atomic oxygen, 
S{0) = 0.5,0.3,0.1 in panels (a), (b), (c), respectively. In all cases, 5(C) = S{N) = 0.1 and S = 1.0 for all other 
neutral species; initially, nc/no = 0.97 w 1 in the gas phase, and n(N)/n(N2) = 18. 



